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Abstract 

The set of world lines for the non-relativistic quartic oscillator satisfying Newtons equation of motion for all space and 
time in 1-1 dimensions with no constraints other than the spring restoring force is shown to be equivalent (1-1-onto) to 
the corresponding set for the harmonic oscillator. This is established via an energy preserving invertible linearization 
map which consists of an explicit nonlinear algebraic deformation of coordinates and a nonlinear deformation of time 
coordinates involving a quadrature. In the context stated, the map also explicitly solves Newtons equation for the 
quartic oscillator for arbitrary initial data on the real line. This map is extended to all attractive potentials given by 
even powers of the space coordinate. It thus provides classes of new solutions to the initial value problem for all these 
potentials. 
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Parti 



INTRODUCTION 

The set of world lines for the non-relativistic quartic oscillator satisfying Newtons equation of motion for all space and 
time in 1-1 dimensions with no constraints other than the spring restoring force is shown to be equivalent (1-1-onto) to 
the corresponding set for the harmonic oscillator. This is established via an energy preserving invertible linearization 
map which consists of an explicit nonlinear algebraic deformation of coordinates and a nonlinear deformation of time 
coordinates involving a quadrature. In the context stated, this result also explicitly solves Newtons equation for the 
quartic oscillator for arbitrary initial data on the real line. No approximations are involved! 

Specifically, each world line for the non-relativistic harmonic oscillator satisfying Newtons equation of motion for all 
space and time in 1-1 dimensions with no constraints other than the spring restoring force satisfies 



d 2 

m dl 2 ~ X ^ = X 



where as usual m=mass, ki = spring constant, and kj_ = m(0 and is given by 



j2E 

x (*) = \ 1~ cos a {t-tmax), (1-2) 
V *2 



for -oo < i < °° and E equals the total energy. 

With the invertible map given in Part III, we shall pair each member of the set of all world lines for the harmonic 
oscillator with one for the quartic oscillator with the same mass and total energy obeying Newton's equation of motion: 
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d 2 

m -2 y (t) = -h y 3 (0- (1-3) 
dt 



The invertible map matches potential energies and momenta, hence energies. 



It is important to comment that an extremal, i.e., each solution of the associated variational problem connecting the 
spacetime events (x a ,t a ) and (xb,ib) for the harmonic oscillator ( (y a ,t a ) and (yb,tb) for the quartic oscillator) is a 
segment of the world line with the same energy. Hence, a fundamental role of our energy preserving linearization map 
is to map the set of extremals in spacetime connecting the initial point (x a ,t a ) and the final point (xb,ib) f° r the linear 
(harmonic) oscillator (ho), 1-1 onto the set of extremals in spacetime connecting the initial point (y a ,t a ) and the final 
point (yb,tb) for the quartic oscillator (qo). 

R. C. Santos, J. Santos and J.A.S. Lima first demonstrated the possibility of linearization of the nonrelativisitic quar- 
tic oscillator to that of the nonrelativistic harmonic oscillator [1]. In [1], they approach the problem of dealing with 
nonlinear "oscillators" via the classical Hamilton-Jacobi equation in contrast to our approach via Newton's equations 
of motion. The most important difference is we have to add a fundamental bookkeeping system to their work; namely, 
the signs of what we call x and y over a cycle have to be matched as well as requiring the harmonic oscillator time i 
and the quartic oscillator time t to progress together in a positive manner. Further, we match momenta which coupled 
with matching the potential energies allows us to match total energies. (In [1], they match potential energies.) This 
is critical since we not only map solutions of the quartic oscillator onto a solution of the harmonic oscillator, but this 
allows us to specify a 1-1 map of extremals in spacetime onto extremals in spacetime for the two systems. 

In Part II it is shown that the world lines of any two harmonic oscillators can be mapped 1-1 onto each other. So we 
can select any harmonic oscillator as representative in this context. The essentials of the map we shall use in Part III to 
linearize the quartic oscillator are illustrated in the development in Part II where for this case both the space and time 
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parts are algebraic. 



In Part III, the set of world lines hence the spacetime extremals connecting the spacetime events (y a , t a ) and tb) for 
the quartic oscillator particle system is shown to be equivalent to the set for the linear (harmonic) oscillator connecting 
the spacetime events (x a , i a ) and (x/,, 4) under an energy preserving map which is a nonlinear algebraic deformation 
of harmonic oscillator space coordinates and a nonlinear deformation of the harmonic oscillator time coordinates. The 
latter is given by quadrature and we shall deal with this point in Part III. Part III ends with the ratio of the period of the 
quartic oscillator (T qo ) to that of the harmonic oscillator period (T/,„) which is shown to be inversely proportional to the 
fourth power of energy. 

In Part IV, a summary of the extension of these results to the hierarchy of attractive potentials given by even powers of 
the space coordinate is given. With respect to [1], the remarks made in Part I generalize to potentials considered in Part 



In Part V, notable work in 1 + 1-Dim is briefly described and higher dimensional applications are referenced. 



Finally in Part VI Concluding Remarks, three other possible applications are briefly described. 



IV. 



Part II 



EXTREMAL MAPPING 




The map to a second harmonic oscillator with mass m and space coordinate x is stated in two parts: 
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(A) 




(2.1) 



This implements the physical requirement that \k2X 2 {t) = \~k.2X 2 (t) matching the potential energies at the two different 
times, coupled with matching of the signs of the space coordinates. 

(B) 

di _ jki 



and 




(2.2) 



This follows by requiring 



dx{i) /di = dx(t)/dt. (2.3) 

Given the matching of the potential energies, the matching of the velocities and the masses of the oscillators for all 
values of Ic2 and £2 implies physically matching the momentum and the kinetic energies at the two different times, i.e. 
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Pho(t) = Pho(t)> E ho = E ho = E. 



Equation (2.2) integrates to 



t tmax — \ ~ if tmax) — - if tmax) (2-4) 
V h 03 



Thus 



x{t) =X max COSC0(i-t max ) <^> x(t) =XmaxCOSG)(t-i max ) (2.5) 



or this map implements the correspondence between Newton's equations of motion 



d 2 d 2 
m-pjx{i) — -k 2 x(i) <^ m-p~x(f) — -k 2 x(t). (2.6) 



Finally, all the world lines of the hat system correspond to all the world lines of the tilde system, so we are free to pick 
any harmonic oscillator as a representative. We designate the hat system as our choice. 
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Part III 



EXTREMAL MAPPING ho^qo 



The map to the quartic oscillator with mass m and space coordinate y is stated in two parts: 



(A) 



'=(2k 2 /k 4 )^x/(x 2 ) l l* 



or 



x=(k 4 /2k 2 ) 1 / 2 (y 2 ) l / 2 y, (3.1) 



where y is the space coordinate of the quartic oscillator and we have used the representation sgn(x) = xl{x 2 ) 1 ! 2 and 
similarly for sgn(y). This implements the physical requirement that \k 2 x 2 {f) = \k 4 y 4 (t) i.e. matching the potential 
energies at the two different times, coupled with matching of the signs of the space coordinates. 



(B) 



| = l/2(2k 2 /k 4 ) l /\x 2 (i))-V\ 



and 

di 



dt (2fe 4 A 2 ) 1/2 (y 2 (0) 1/2 , (3.2) 



which results by requiring 
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dx{i)/di = dy{t)/dt. 



(3.3) 



Given the matching of the potential energies, the matching of the velocities and the masses of the oscillators for 
all values of kj_ and implies physically matching the momentum and the kinetic energies at the two different 
times, i.e. pi w (i) = p qo (t), Eho = E qo = E. It is interesting to note that (2.2) is also equivalent to requiring dtldi = 
[(— fcty 3 (f))/(— k2x(t))]~ l = the inverse of the ratio of the forces. 

This map implements the correspondence between Newton's equations of motion 



m-^x{i) = -k 2 x(t) m-j^y(t) = -k4y 3 (t) (3.4) 

in the following straight-forward way. Differentiation of (2.1) shows that (2.1) and (2.2) are self-consistent. Differen- 
tiating (2.1) a second time, invoking (2.2) and assuming the validity of one side of (2.6) then yields the validity of the 
other side of the equivalence in (2.6). 

If one plots the /zo-potential and the go-potentials vs space coordinates, then the horizontal lines in such a plot allow 
one to graphically read off the corresponding coordinates. Because of reflection symmetry in the vertical axis, one 
need only plot the potential for the positive space coordinate. Further, for some y-coordinate near enough to zero but 
not zero, the go-potential is less than the /zo-potential for all nonzero ^-coordinates less than that value. This feature is 
illustrated in Figure 1. 
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Figure 1 : Comparison of Potentials 



In [1], they approach the problem of dealing with nonlinear "oscillators" via the classical Hamilton-Jacobi equation in 
contrast to our approach via Newton's equations of motion. Another way that our work differs from theirs is that we do 
not identify the solutions of the HJ equation, rather we only assume their existence. Another difference is we have to 
add a fundamental bookkeeping addition to their work; namely, the signs of what we call x and y over a cycle have to 
be matched as well as requiring the harmonic oscillator time t and the quartic oscillator time t to progress together in a 
positive manner. Further, we match momenta which, coupled with matching the potential energies, allows us to match 
total energies. (In [1], they match potential energies.) This is critical since we not only map solutions of the quartic 
oscillator onto a solution of the harmonic oscillator, this allows us to specify a 1-1 map of world lines in spacetime 
onto world lines in spacetime for the two systems. 



(y(t),P go (t)) 




Figure 2: Phase Space 



The corresponding phase space diagrams are interesting, but very case dependent in terms of and £/,„. The hori- 
zontal coordinate scales are set by x max = \f2E~jk~2 and _y max = \J^E jk/\. As a result y max can be less than, greater than 
or equal to x max . The vertical scales are set by the equality of the linear momentum at the different times, i.e. m c -^(t) = 
m4fX(t) which also implies that = m^x . If one were to draw a phase diagram, horizontal lines yield the 

max max 

corresponding phase space coordinates (x(f), %(f)} -O- (y(t), ^ (0)- Thus, while both make one complete cycle together 
in terms of phase coordinates, they do it in different times. This point is important to emphasize; one complete cycle 
for the ho in the phase space diagram corresponds to one complete cycle for the qo in the phase space diagram, but this 
is done in general in different periods. These features are illustrated in Figure 2 fory max < x max . Because of symmetry, 
one need only graph the first quadrant of the phase diagram. 



Now, the nonlinear deformation of time given by (3.2) is in terms of a definite integral. 



(1/2) (2k 2 /k 4 ) 1 / 4 f 



'b 1 



(*2(f))V4 



dt = tj, — t a 



(3.5) 
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Equation (3.5) is the operative quadrature that establishes the connection between the two "times" t and i . Mathematica 
[2] provides a complete program for evaluating this integral. 

Now we turn to evaluating (3.5) for the period of the quartic oscillator x qo in terms of the period T/, of the harmonic 
oscillator. It is sufficient and convenient to set t max = n/2 in (1.2) or take the world line corresponding to 

x(i) = x max sincof, (3.6) 

where as before x max = \J (2E/k2). Then (3.5) becomes for one cycle of the quartic harmonic oscillator (or linear 
oscillator) 

1/2 {k\lUE) l ' A r° ^ dt = x qo . (3.7) 

J ° ( (sin 2 ©f)i/4j 

Equation (3.7) is an elliptic function of the first kind [2]. 

A simplification is possible, namely the change of variable 9 = at. Equation (3.7) then becomes 



\ V 4 f 2K 1 

1/2 [l4/k 4 E] a' 1 I - f ^ dd = T qo , (3.8) 

(sin 2 ©) 1 / 4 



7 jo 



where, again, o = 27zVT/ ! „. 
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Equation (3.8), using symmetry, yields the period T qo of the oscillator in terms of the period T), as follows 




(3.9) 



where using Mathematica [2] the term in square brackets equals .83. It follows from (3.8) that for fixed k^ and kt, we 
have X q(> IXho ~ £~ 1/4 



In this section we present a summary of the extension of these results to the hierarchy of attractive potentials given by 
even powers of the space coordinate paralleling that given in Part II - Part III. 

First, we shall outline the mapping of the harmonic oscillator extremals onto the extremals of a each member of an 
hierarchy of attractive oscillators with coordinates y2 n (t)\ n = 2, 3, 4, ... characterized by even positive power law 
potentials. (The case, y = yA, which is included in the hierarchy, has been the subject of the preceding paragraph.) In a 
straight forward manner the mappings in Part III, generalize and yield the following relationships: 



Part IV 



EXTREMAL MAPPING 




HIERARCHY 
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(A) 



y 2 „ = (nfe2A2„) 1/2 "^ 2 ) (1/2)(1 - fl)/ " 

x={k 2n /nk 2 ) l l 2 y 2n {yl n )^l\ (4.1) 



which is the generalization of (3.1). The generalization of (3.2) is given by: 



(B) 



% = n-^- i y 2n {k 2 /k 2n ) l i 2n (x 2 (i))- {n - l) i 2n (i), 

at 



and 



f t =Vn(k 2n /k 2 )^(yl n )^y\ (4.2) 



These mappings take the space-time extremals of the linear oscillator with coordinates (x, t) and map them onto the 
space-time extremals of the 2n th oscillator with coordinates (y 2n , t). 



A straightforward calculation yields 



d 2 d 2 
m-j2jx(i) = -k 2 x(i & m-^y 2n (t) = ~k 2n yf n l (t). (4.3) 



2n-l, 



Further, as a consequence of the above, we have the following equality for the conserved total energies 
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The operative deformation of time given by (4.2) becomes in integral form 



(4.4) 



t b -t a = [ b n- {2n -^l 2n {k2/k 2n ) l l 2n { X 2 {t))- {n -^l 2n di (4.5) 
This is implemented using Mathematica [2]. 

All of the analyses presented in Part III can then be extended to the members of the hierarchy. 

PartV 

NOTABLE CONSERVATIVE 1 + 1-DIM WORK 

Whittaker [3, p. 64] solved all conservative 1 + 1-Dim mechanical systems up to a quadrature by solving the energy 
relationship for the velocity and integrating: 

t b -ta= r±V2m^E -V(y')dy'. (5.1) 

ha 

For the quartic oscillator, we can cast (5.1) into the following form 
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2jm&M)& de ' (5 - 2) 



where 



(y 2 ) l / 2 y = sinO. (5.3) 
We have accounted for the sign of the velocity with cosd. 

Similar forms can be derived for the rest of the attractive hierarchy studied in this paper. 

One of the most interesting formulations is the one found by a "regularization" scheme for the differential equations 
[4, p. 14-17] describing the negative energy 1 + 1-Dim Kepler problem. Interestingly, this leads to a cycloid solution 
and related harmonic oscillator solution. This problem has an old and distinguished past and the author has not done 
justice to it here. This solution is not obtainable by matching of potential and kinetic energies as employed in this 
paper. The utility of the 1 + 2-Dim and the 1 + 3-Dim versions of their regularization in the Feynman Path Integral 
Method is prominent in the work of [5], [6]. 

The radical solutions of 1 + 2-Dim central force problems with angular momentum as one constant of the motion and 
energy as the other are addressed by Whittaker [3, p. 80-81]. He solves them in the spirit of the 1 + 1-Dim problem 
cited above, but in two dimensions it requires two quadratures. 




15 



Part VI 



CONCLUDING REMARKS 

Physically what we have done with our linearization map is to view the tape of the space-time evolution of these non- 
linear systems using the right optical lens to accomplish the nonlinear deformation of space and playing each frame 
according to a nonlinear deformation of time. When we do this we see that the quartic oscillator evolves like an har- 
monic oscillator. 

The mappings given in Parts III-IV provide new classes of exact solutions for nonlinear spring systems. 

The essence of the maps presented here is that they involved two conservative mechanical systems and if it makes sense 
they match the potential energies and momenta of the two systems. Within this context, its utility arises if solutions (a 
solution) of one of the systems are (is) known. This means e.g., that one could linearize the hierarchy 

or one could treat each motion of a quartic oscillator bouncing off an infinite barrier as a sequence of maps to a harmonic 
oscillator capturing the discontinuous velocity at the barrier, or one could map the radial equation in the procession 
of Mercury onto that of the inverse law force (4), [7, p. 194]. The latter involves extracting the zeros of polynomials 
cubic in ji = j for the algebraic part of the map. This is a nonlinear problem to nonlinear problem in the context of 
boundstates. 

Special thanks go to my Department of Mathematics colleague Robert Varley. He spent enumerable hours over a four 
year period of time discussing this work with me. His comments, questions and posing of challenging related problems 
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helped to clarify for me many aspects of this work. 

The author wishes to aknowledge the referee's role in particular, in directing him to the work in [5] and [6]. 

The author wishes to thank Professor Howard Lee for insightful discussions and his constant encouragement. The idea 
to emphasize the quartic oscillator was his. 
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